Critical spin transport in Bose gases 
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We consider spin transport in a two-component atomic Bose gas in three dimensions, at tem- 
peratures just above the critical temperature for Bose-Einstein condensation. In these systems the 
spin conductivity is determined by spin drag, i.e., frictional drag between the two spin components 
due to interactions. We find that in the critical region the temperature dependence of the spin 
conductivity deviates qualitatively from the Boltzmann result and is fully determined by the critical 
exponents of the phase transition. We discuss the size of the critical region where these results may 
be observed experimentally. 
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Introduction. — The research field of spin electronics 
or spintronics, concerned with practical applications of 
the electron spin, has renewed interest in spin currents 
[T]. In part as a result of these efforts, it is now un- 
derstood that there are several fundamental differences 
between charge currents and currents of spin angular mo- 
mentum. For example, the latter are even under time 
reversal-symmetry operations and can thus in principle 
flow without dissipation in ordinary conductors [H [3], 
contrary to electric currents. Furthermore, the charge 
conductivity is infinitely large in Galilean invariant sys- 
tems, whereas spin currents can then still decay due to 
spin-drag effects, i.e., friction between two spin states 
due to interactions [H |S] . Finally, spin and charge cur- 
rent couple in a completely different way to other de- 
grees of freedom in the system, most notably order pa- 
rameters such as the magnetization or a superconduct- 
ing condensate. For example, a spin current can exert a 
so-called spin transfer torque on the magnetization of a 
ferromagnet [5HS] , a phenomenon that is currently inten- 
sively studied in part because of its promise for magnetic- 
memory applications. On a more fundamental level, 
there have been several studies on the interplay between 
spin currents and the critical fluctuations in the magne- 
tization that occur close to the Curie temperature for the 
ferromagnetic phase transition [lOj . Such magnetic phase 
transitions form, together with superconducting phase 
transitions, the overwhelming majority of phase transi- 
tions occurring in electronic condensed-matter physics. 

In this Letter we consider the effect of another phase 
transition, i.e., Bose-Einstein condensation, on spin 
transport. The system we consider is a spin mixture of 
trapped ultracold bosonic alkali atoms that differs in sev- 
eral important ways from electronic solid-state systems. 
First, the particles are bosons rather then fermions (elec- 
trons). Second, cold- atom systems are disorder free and 
hence the only contribution to the spin conductivity is 
the above-mentioned spin-drag effect. The atomic in- 
teractions which are responsible for this drag are short 
ranged as opposed to the Coulomb interactions between 
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FIG. 1: (Color online) Spin-drag relaxation rate l/rsd at con- 
stant temperature as a function of distance from the critical 
point expressed in terms of the chemical potential difference 
/i — ^c- Lines 1, 2 and 3 represent the relaxation rate for 
a/Ath = 9 X 10"^, 2(9 X 10"^) and 3(9 x 10"^), respectively, 
where a is the scattering length and Ath = \f2TvW]JrnkBT) is 
thermal de Broglie wavelength. Upon approaching the tran- 
sition from above, the spin-drag relaxation rate shows an up- 
turn due to Bose enhancement that is ultimately completely 
suppressed by fluctuations in the critical region. The dotted 
line 4 represents the Boltzmann result that does not include 
critical fluctuations [T^. The small quantitative difference 
between the lines 1-3 and the Boltzmann results far from 
criticality arises because we neglect vertex corrections in the 
calculations that lead to the curves 1-3. 



the electrons. 

Spin drag in bosonic cold-atom mixtures was recently 
studied by two of us using an approach based on the 
Boltzmann equation. It was found that the bosonic na- 
ture of the particles lead to an enhancement of spin-drag 
effects at low temperatures [11]. This should be con- 
trasted with Fermi-liquid behavior that as a result of 




FIG. 2; Sunset diagram for the self-energy of spin-up particles 
(the diagram for spin-down is obtained by reversing all spins). 



Pauli blocking leads to suppression of interaction effects 
at low temperatures so that the so-called spin-drag relax- 
ation rate l/TgdiT), which is equal to the inverse of the 
spin-transport relaxation time Tsd(T), vanishes quadrat- 
ically with temperature T |121 113) except in the vicinity 
of a superconducting |14j or ferromagnetic phase transi- 
tion '15]. For bosons close to the critical temperature 
for Bose-Einstein condensation, it was found that the 
Boltzmann approach incorporates the phase transition 
at the mean-field level and gives l/rsd(T) — l/Tsd(T'c) ~ 
-1/C(T) - Tc - r J^, where S.(T) is the correlation 
length that diverges at the phase transition. 

Our main findings are presented in Fig. IT] which shows 
the spin-drag relaxation rate as a function of the dis- 
tance to the critical point, determined by the difference 
of the chemical potential fi from its critical value ^c- 
The dotted line shows the Boltzmann result (described 
above) and the solid lines are the results found using an 
approach based on the Kubo formula and approximating 
the atomic self-energy with the so-called sunset Feynman 
diagram shown in Fig.|2] The Hartree diagram should in 
principle also be included in the self-energy but since this 
can be achieved by a simple redefinition of the chemical 
potential we do not consider it here. Within the lat- 
ter approximation we find that the spin-drag relaxation 
rate qualitatively agrees with the Boltzmann result for 
temperatures not too close to the critical temperature, 
but in the critical region deviates and goes to zero at 
the phase transition according to l/Tsd{T) ^ 1/£,{T). As 
we discuss in detail below, an exact scaling ansatz con- 
firms that the spin-drag relaxation rate vanishes, and in 
terms of the critical exponents z, 77, and v we find that 
^/Tsd{T) ^ ii^z-d+2-2ri^ Below wc also discuss the size 
of the critical region where these effects can be measured. 

Spin-drag conductivity. — We consider a three- 
dimensional homogeneous gas of bosonic atoms of mass 
m, with two spin states that couple with opposite sign 
to an external force F. As we discuss in more detail 



below, this force can in a cold-atom experiment be im- 
plemented by a magnetic-field gradient. This force leads 
to a nonzero spin current js according to js — crgF, 
where a^ — uTsd/m is the spin conductivity in terms 
of the spin-drag relaxation time and the density n per 
spin state (we consider here only the balanced case 
where the densities of the two spin states are equal). 
The Kubo formula for the spin conductivity is given 
by (Ts = — lim[j_>.o 3tn Hiii, {k = 0,uj) juj in terms of 

the Fourier transform of the retarded spin-current spin- 
current correlation function n^y(a;, t; a;', t') = i9{t — 
t'){[jj,'{x,t),j^{x',t')])/h, where the expectation value 
(• • •) is taken in equilibrium and 
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In the above expression the bosonic Heisenberg annihi- 
lation operator is denoted by ^{x,t) with a € {t;i} 
labeling the spin states, and we note that a takes the 
respective numerical value -I- or — if it is not used as a 
label. 

To evaluate the correlation function in the Kubo for- 
mula we neglect interactions between atoms of like spins 
that are of minor importance for spin-drag effects. The 
hamiltonian then reads 
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da;?Aj(a;,t)-0|(a;,t)i/>^(a;,i)Vi;(a;,i) , (1) 



with r^^ = Anh^a/m the two-body T matrix in terms of 
the interatomic scattering length a between the two spin 
states. 

Upon ignoring vertex corrections to the correlation 
function II},^ (fc,w) the Kubo formula for the spin con- 
ductivity is worked out to yield 
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{k,Lj) /nh the spectral 
function for particles with spin a, obtained from 
the Fourier transform of their retarded Green's func- 
tion that is defined by Ga {x,t;x',t') = i6{t — 

t')( 'il}a{x,t),i}\{x',t')\. Furthermore, Nsihuj) = 
l^gPHu} _ 2J-1 jg ^]-^g Bose-Einstein distribution function, 
with /3 — l/fceT the inverse thermal energy. 

The lowest-order (in interatomic interactions) diagram 
that gives a finite conductivity is the sunset diagram in 
Fig. [2] for the atomic self-energy 

Jm[hj:i+Hk,Lj)] = 



^(rr.2B^2 f dk2 f dk3 



(27r)3 



3500 



X (27r)^(5(fc + ^2 - fca - k4)5{uj + fi + ek^- ^h ~ ek^) 

x[NB{k2){l + iVs(fc3))(l + Neiki)) 

-{l + NB{k2))NB{k3)NB{ki)], (3) 

with Efc = fv^k"^ /2m. The imaginary part of the self- 
energy is calculated numerically from Eq. (pi). Its real 
part is obtained by using a Kramcrs-Kronig relation. 
The spectral function then follows by using Ga {k^uj) — 
h/{fibj ~ ek + fJ- — hY^a {k,iij)). This real part shifts the 
critical chemical potential from its noninteracting value 
of zero to the positive value fic — 9^e [/iS^+-'(0, 0)] . 

The result for the spectral function is shown in Fig. [3] 
as a function of frequency. For a given momentum, the 
spectral function exhibits a rather sharp Lorentzian peak 
corresponding to a quasi-particle excitation. For such 
a Lorenzian spectral function, the frequency integral in 
Eq. ([2| can be performed and the resulting conductivity 
is then found to be proportional to the life-time t(A;) = 
— S/(23m[fiE(+)(fc,a;fc)]) of the quasi-particle, where Wfc 
is the solution of Huk = e^ — ^ + 9^e[fiS(+)(fc, Wfc)]. The 
evaluation of the expression for the spin conductivity 
in Eq. ([2]) with the above expression for the self-energy 
leads to the results shown in Fig. [l] for various scattering 
lengths. 

To gain more insight in the breakdown of the Boltz- 
mann approach in the critical region, we consider the 
spectral function at low momentum inside (outside) the 
critical region, corresponding to the left (right) peak in 
Fig. [3] The vertical lines correspond to efc — /i. The peak 
in the spectral function is shifted considerably from its 
non-interacting value ey. — ^ va the critical region. Since 
the Boltzmann approach does not take into account the 
shifts in the quasi-particle energy beyond first order in 
the interaction it does not capture the shift in the crit- 
ical region correctly. The latter is crucial for the diver- 
gence of the conductivity. The importance of the real 
part of the atomic self-energy in the critical region is also 
demonstrated by its importance in determining the up- 
ward shift in the critical temperature due to interactions, 
that is correctly found to be of order O (a/Ath) p^ . 

Critical phenomena and scaling. — From the calcu- 
lation that is based on the sunset diagram for the self- 
energy we find numerically that the spin conductivity 
diverges as Us ^ l/^/ic — /i. This result is understood 
on a more general level by considering scale invariance of 
the system. Near criticality, we have as a result of the 
simple (linearized) renormalization-group flow near the 
fixed point that the spectral function scales as [T7] 
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where A is an arbitrary dimensionless scaling parameter 
and z, V, and rj are critical exponents. The relation be- 



FIG. 3: The spectral function at low momentum inside (left) 
and outside (right) the critical region. The dashed and dotted 
lines show the value of e^ — p for the spectral function inside 
and outside the critical region, respectively. 



tween the correlation length ^ and the chemical potential 
is ^ ^ 1/l/i — ^ic\'^ . From this scaling ansatz and Eq. (pi) 
we find that ct^ - ^^-^+2-2,, -- Ia^ - ^Jic\^''^^^'^+^^'^, 
with d the number of spatial dimensions. 

For the sunset diagram in three dimensions, we have 
v — 1/2, z = 2, and 77 = 0, in agreement with the numer- 
ical results. It is interesting to note that the behavior of 
the spin conductivity and the spin-drag relaxation time 
depends not only on the static critical exponents v and 77 
but also on the dynamical exponent z. We also note that, 
even though we ignored interactions between atoms with 
the same spin in our perturbative calculation based on 
the Feynman diagram in Fig. [2] the results based on the 
scaling ansatz are exact close to the critical temperature 
and do include these interactions. 

Following the reasoning of Hohcnberg and Halperin 
;19] the factor ^^-d ^^ these results is understood as fol- 
lows. A spin-dependent force acting on a region with 
(fluctuating) spin density n^ is balanced by viscous forces 
so that ns£,^F ^ S^^-q^^Vs/S,'^, where Vg is the spin velocity 
and rjy the viscosity. Using that jg — n^Vs = cf^F this 
yields as ^ (.'^ (n'j,) / r]y . We have that (n^) - Xs/^'^ }19j . 
with Xs the spin susceptibility. In order to obtain full 
agreement with our result found from the scaling ansatz 
we thus need to have that the ratio Xs/Vv ^ £,^~'^'^. In 
future work we intend to investigate this conjecture in 
more detail. 

Discussion and conclusions. — We have incorporated 
the effect of critical fluctuations on the behavior of the 
spin-drag relaxation rate near the critical temperature 
for Bose-Einstein condensation. We found that the en- 
hancement of the spin-drag relaxation rate due to Bose 



enhancement of interatomic interactions, predicted by 
the Boltzmann equation, is suppressed by critical fluctu- 
ations sufficiently close to the critical temperature. Nu- 
merically, we found the critical region to be proportional 
to the square of scattering length |A^| « 60(a/Ath)^- 
An estimate based on the Ginzburg criterion [5T] con- 
firms this result. Hence, the size of the critical region 
may be enlarged by increasing interatomic interactions 
near a Feshbach resonance. Furthermore, the Ginzburg 
criterion in d dimensions leads to |A/i| ~ (o/Ath)^^^'*~'''. 
Reducing the dimensionality of the system therefore also 
increases the critical region. With respect to these re- 
marks it is important to note that recent experiments 
with ultracold bosonic atoms have succeeded in access- 
ing the critical region and measuring the exponent v |20| . 

The spin conductivity and spin-drag relaxation rate 
can be measured directly in a drag measurement in which 
the two clouds of different spin feel a different force due a 
magnetic-field gradient. Another method is to study the 
damping of the spin-dipole mode that is fully determined 
by the spin-drag relaxation rate. 

The main approximation leading to our results is to 
neglect vertex corrections in the evaluation of the spin- 
current spin-current response function. The Boltzmann 
equation is known to include vertex corrections that es- 
sentially lead to a replacement of the single-particle relax- 
ation time by the appropriate transport relaxation time. 
In the absence of exact cancelations, which we do not ex- 
pect to occur for the spin-drag conductivity, there is only 
a quantitative difference between these two time scales, 
and we attribute the difference between our Kubo ap- 
proach and the Boltzmann approach sufficiently far away 
from the critical region to be due to this difference in time 
scales and hence to be due to neglecting vertex correc- 
tions. This assumption is strengthened by noting that 
far away from the critical region the Kubo and Boltz- 
mann approach have qualitatively the same temperature 
dependence. 

In this work we have considered the effect of ther- 
mal critical fluctuations since the phase transition to 
the Bose-Einstein-condensed state takes place at nonzero 
temperature. An interesting direction for future work is 
to investigate also the influence of the vicinity of a quan- 
tum critical point |.22j on spin transport in Bose gases, for 
example by considering the system in an optical lattice 
where the Mott-insulator-to-superfluid quantum phase 
transition occurs [23]. An additional interesting feature 
of this system is that due to the presence of the optical 
lattice, which breaks Galilean invariance, now also charge 
(mass) transport can be considered. 
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